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Abstract 

The dominant production mechanism for be bound states in high energy processes 
is the production of a high energy b or c quark, followed by its fragmentation into the 
be state. We calculate the fragmentation functions for the production of the S-wave 
states B c and B* to leading order in the QCD coupling constant. The fragmentation 

probabilities for b — > B c and b —* B* are approximately 2.2 x 10 -4 and 3.1 x 10 , 
while those for c — > B c and c — > B* are smaller by almost two orders of magnitude. 



The study of the charmonium and bottomonium systems have played an important 
role in the development and eventual acceptance of quantum chromodynamics (QCD) as 
the field theory of the strong interactions. Bound states containing top quarks may not 
exist as identifiable states, because the widths of the states due to top quark decay may be 
larger than the splittings between the energy levels. Thus the final frontier for the study 
of heavy quark- ant iquark bound states may be the be system. To assess the prospects for 
the experimental study of be mesons in present and future colliders, it is important to have 
accurate predictions for their production rates. The dominant production mechanism for 
high energy be mesons is the fragmentation of b and c quarks. The fragmentation process 
is described by a universal fragmentation function D(z) that gives the probability for the 
splitting of the parton into the desired hadron plus other partons. For a hadron composed 
of only heavy quarks, this fragmentation function can be calculated using perturbative QCD 
U|. The fragmentation functions for the splitting of gluons and charm quarks into S-wave 
charmonium states have been calculated to leading order in the QCD coupling constant a s 
[p], 0. In this paper, we calculate the fragmentation functions for the splitting of b and c into 
the l So ground state B c and into the 3 Si state B*. These universal fragmentation functions 
can be used to calculate the direct production rate of B c and B* in any high energy process. 

We begin by summarizing previous calculations of the production rates for be mesons. 
A perturbative QCD calculation of the cross section for e + e~ — > B*bc was carried out by 
Clavelli in 1982 |J. The cross sections for e + e~ — ► B c bc and e + e" — > B* be were subsequently 
calculated incorrectly by Amiri and Ji |4j] and then correctly by Chang and Chen ||. Chang 
and Chen followed Amiri and Ji in expressing their results in terms of fragmentation functions 
D(z), but they did not realize that these fragmentation functions are universal and can also 
be applied to the hadronic production of B c and B*. The total production rates for B c in 
e + e~ and hadron colliders, including those that are produced by cascade decays from B* and 
other bc bound states, have been estimated using a modified version of the parton shower 
Monte Carlo HERWIG |J. The modified version included two new parameters describing 
the probabilities for the production of cc pairs, and the results were sensitive to the values of 
the bottom and charm quark masses. The ratio of the cross section for B c production and b- 
quark production was estimated to be 1.0 x 10~ 3 for e + e~ collisions at LEP and 1.3 x 10 -3 for 
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pp collisions at the Tevatron, consistent with a previous crude estimate |?j of the probability 
for a b to form a be meson. Given the adjustable parameters in the Monte Carlo calculation 
of Ref. 0, it is difficult to estimate the errors in its predictions. It is clear that more reliable 
calculations of the production rates of be mesons would be valuable in assessing the prospects 
for the discovery and investigation of these new particles. 

The fragmentation contribution to the differential cross section for the direct produc- 
tion of a B c of energy E in any high energy process is the term that survives in the limit 
E/rrih — > oo. It has the general form 

da(B c (E)) = t dz da(t(E/z),fi) D^ Bc (z, ^) , (1) 

where the sum is over partons of type % and z is the longitudinal momentum fraction of the 
B c relative to the parton. The physical interpretation of (|J) is that a B c of energy E can 
be produced by first producing a parton % of larger energy E/z which subsequently splits 
into a B c carrying a fraction z of the parton energy. The expression (|I|) for the differential 
cross section has a factored form: all the dependence on the energy E is in the parton 
subprocess cross section a, while all the dependence on the heavy quark masses m c and m;, 
is in the fragmentation function Di_>B c - The dependence on the arbitrary factorization scale 
// cancels between the two factors. Large logarithms of Ej\i in the subprocess cross section a 
can be avoided by choosing \i on the order of E. Large logarithms of ju/mj then necessarily 
appear in the fragmentation functions Di^ Bc (z, /x), but they can be summed up by solving 
the evolution equations || 

^-^-Di^ Bc (z,ii) = — P*-+j(z/y,n) Dj->B e (y,li) , (2) 

ofi . j z y 

where -Pi->j (x, //) is the Altarelli-Parisi function for the splitting of the parton of type % into 

a parton of type j with longitudinal momentum fraction x. For example, the b ^b splitting 
function for a bottom antiquark with energy much greater than its mass is the usual splitting 
function for quarks and antiquarks: 

= ^(tzt) + > ( 3 ) 
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where f(x)+ = f(x) — 5(l—x) Jq f(x')dx'. The boundary condition on the evolution equation 
(0) is the initial fragmentation function Di^B c {z, /io) at some scale /io of order m b . The initial 
fragmentation function can be calculated perturbatively as a series in a s (/io). At leading 
order in a s , the initial fragmentation function for b — > B c can be calculated using the 
expression 

D- b _ Bc (z) = -— ds6\s hm , (4) 

167r 2 7 \ 2 1 — z) 9o-*t» |yvior 

where .M is the matrix element for the production of a B c and c with total 4-momentum q 
and invariant mass s = q 2 and .Mo is the matrix element for the production of a b with the 
same 3- momentum q. If the momentum of the B c is p = (po,Pi;P2,P3) in a frame where q = 
(go, 0, 0, qs), then the longitudinal momentum fraction z is defined by z = (po + P3) / {qo + Q3) ■ 
The limit go — > 00 in (§) can be taken in any Lorentz frame where ga/go — > 1. 

We proceed to calculate the matrix element yVl for the production of B c in the frag- 
mentation limit go >> m b . If g = (g , 0, 0, g 3 ), it is convenient to carry out the calculation 
in the axial gauge associated with the 4- vector n = (1, 0, 0, —1). In this gauge, the fragmen- 
tation contribution comes only from the Feynman diagram for b* — > bcc shown in Figure 1. 
In covariant gauges, there are also contributions from other diagrams but Ward identities 
can be used to manipulate the sum of the diagrams into an expression equivalent to that 
obtained directly from Figure 1 in axial gauge. The amplitude for Figure 1 in this gauge can 
be reduced to 

= g^EEL- 1 2x? f(2U+m b ) (/-2M) 75 
3v / 3TMr(s -mf) 2 V 1 KF 1 

' ' (p* + M) l5 yl) v(p>) , (5) 



n ■ (q — (1 — r)p) 



where r = m c j (m b + m c ), M = m b + m c , and R(0) is the nonrelativistic radial wavefunction 
at the origin for the B c . The Dirac spinor f is the matrix element for the production of a b of 
momentum q = p+p f . Its explicit form is not needed to calculate the fragmentation function 
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). Summing over final spins and colors, the square of the matrix element (0) reduces to 



647rct 2 | J R(0)| 2 1 



where A is a Dirac matrix that depends on q and p. We need only keep the terms in A 
that are of order m^q . While q and p both have components of order q , the invariant 
mass s = q 2 is of order m 2 in the fragmentation region. Simplifying the Dirac matrix A by 
dropping terms which are suppressed by powers of m^/qo, it reduces to 

A = (s 2 - 2(2 - 4r + r 2 )M 2 s + (1 - r)(3 - r)(l - 4r + r 2 )M 4 ) i 

- (s - (2 - r) 2 M 2 ) (s - (1 - r) 2 M 2 ) i> 

+ 2 ( S -(l-r) 2 M 2 ) 2 (^) 

n ■ q — (1 — r)n ■ p 

~ 2 n . q -' l { Zr)n- p M2(s ~ (1 " TfM ' ] ((1 " 2r) 4 ~ (1 " r) A 

+ -^-h-rt (,-(!_, fMT/. (7) 
(n - g — (1 — rjn ■ p) 2 

Since the coefficients of $ and ^ in (^) are all manifestly of order m 4 , we can make the 
substitution p = zq which is accurate to leading order in m b /q . The Dirac trace in (0) 
is then proportional to tv(TTjf). The square of the matrix element for the production of 
a real b of momentum q is |.Mo| 2 — tr(IT(j/ — m^)), which reduces in the fragmentation 
limit to tr(rr^). Dividing \Ai\ 2 in ([]) by |.Mo| 2 and inserting it into (|j), we obtain the 
fragmentation function as an integral over s: 



*_>B C V ; " 277rm 3 7 ^ X _ z j ^( 1 _( 1 _ r ) /2 )2( s _ m 2)2 

(2(1 - 2r) - (3 - 4r + 4r 2 )^ + (1 - r)(l - 2r)z 2 )rM 4 4r 2 (l - r)M 6 \ 
(1 - (1 -r)z){s-m 2 Y (s - m 2 ) 4 J ' 



(8) 



Evaluating the integral over s, we obtain our final expression for the fragmentation function 
at the initial scale /xq: 

2a s (u ) 2 | J R(0)| 2 rz(l-z) 2 ( , 9n 9 

D^ ^/io = ^° 3 ?! 71 ^ 6 - 18 l-2r z + 21-74r + 68r 2 z 2 

8l7rm2 (1 — (1 — r)z) b V 



4 



- 2(1 - r)(6 - 19r + 18r 2 )z 3 + 3(1 - r) 2 (l - 2r + 2r 2 )z 4 J , (9) 

where r = m c / (m&+m c ). This agrees with the result of Chang and Chen ||] up to the overall 
normalization. Chang and Chen normalized their fragmentation functions so that they 
integrated to 1: Jq dzD(z) = 1. Setting m& = m c in (Q), we also recover the fragmentation 
function for c — > rj c calculated in Ref. Q. The incorrect result of Amiri and Ji [|J can be 
obtained by keeping only the term proportional to l/(s — m 2 ) 2 in the integrand of (j8|). 

The initial scale /j,q in the fragmentation function and in the running coupling constant 
in ([|) should be chosen on the order of and m c . We would like a more specific prescription 
for fio that remains physically sensible even when m c / nib is much less than or much greater 
than 1. In the Feynman diagram in Figure 1, the invariant mass of the gluon is r(s — ml) 
and the virtual b is off shell by s — m 2 . For generic values of z, s is of order m 2 , and the b 
and the gluon are off-shell by amounts of order (m& + m c ) 2 and m c (mb + m c ), respectively 
However a large part of the fragmentation probability comes from the region where s is 
near its minimum value (m& + 2m c ) 2 and the virtual b and the gluon are only off-shell by 
4m c (mb + m c ) and 4m 2 . As a compromise between the various scales in the problem, we 

take fj,Q = \J 4m c (mb + m c ) . In Ref. ||, Chang and Chen used the value a s = 0.15 for the 

running coupling constant, which corresponds to a scale fiQ of order Mz- As shown by our 
analysis, the scale for the fragmentation process is definitely set by the heavy quark masses. 
Thus Chang and Chen underestimated the rates for production of B c and B* in Z° decay 
by about a factor of 2. 

The fragmentation function for a b to split into the 3 Si state B* can be calculated in 
the same way as for B c . The matrix element for the production of B*c is 

3V37TM r(s - ml) 1 \ 

Q 1 11 ^ \ 

+ — ( — tt^t U+M)-frL) «(p0 . (io) 

n • [q — (1 — r)p) J 

where e M (p) is the polarization vector for B*. Following the same path as in the B c calcula- 
tion, we obtain an expression for the fragmentation function as an integral over the invariant 
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mass s: 



4a 2 s \R{0)\ 2 f , n f M 2 m 2 \ /(l - z)(l + 2rz + {2 + r 2 )z 2 )rM 2 
D- b ^ B Jz) = —— — ds6\s 



27nm^ J \ z 1 — z I I (1 — (1 — r)z) 2 (s — m 2 ) 2 
(2(1 + 2r) - (1 + 12r - Ar 2 )z - (1 - r)(l + 2r)z 2 )rM 4 12r 2 (l - r)M QS 



(1 - (1 - r)z)(s - ml) 3 (s~m 2 ) 4 , 

Integrating over s, the final result for the fragmentation function for b to split into B* is 



.(11} 



2o; s (uo) 2 | J R(0)| 2 rz(l-z) 2 , , , 

D- b ^(z^ ) = '^JJ n (1 _ (1 _1J 2 ~ 2(3 -2r)z + 3(3 - 2r + 4r 2 )z 2 



c 



- 2(l-r)(A-r + 2r 2 )z 6 + (1 - r) 2 (3 - 2r + 2r 2 ) z* J . (12) 

This agrees with the result of Chang and Chen [0 up to the overall normalization. Setting 
vn\j = m c in QT2]), we recover the fragmentation function for c — > ip calculated in Ref. H. The 



incorrect result of Amiri and Ji M can be obtained by keeping only the term proportional 
to l/(s — ml) 2 in the integrand of flTTl) . 



The fragmentation functions for a charm quark to split into B c and B* can be obtained 
from (^) and (|T2"|) by interchanging the masses m^ and m c : 

, /N 2ct s (u') 2 |i?(0)| 2 (l-r)^(l-^) 2 / , 9x „ 

LUb^X = 3 n v> 6 + 18(1 -2r> + 15 - 62r + 68r 2 z 2 

8l7rmg (1 — rzy \ 

- 2r(5 - 17r + 18r 2 )z 3 + 3r 2 (l - 2r + 2r 2 )/ ] , (13) 



- 2r(5 - 3r + 2r 2 )^ 3 + r 2 (3 - 2r + 2r 2 )^ 4 ) , (14) 



where r = m c / (m^ + m c ) and fi' = ^jAm^mi, + m c ). 

The fragmentation functions for b — > B c and 6 — > B* are shown in Figure 2. The 
radial wavefunction at the origin for the B c and B* has been estimated from potential 
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models to be |i?(0)| 2 = (1.18 GeV) 3 JT0|. For the quark masses, we use m& = 4.9 GeV and 
m c = 1.5 GeV. The solid curves in Figure 2 are the initial fragmentation functions at the 
scale /J = 6.2 GeV, where the running coupling constant has the value a s (fi) = 0.20. The 
dotted curves are the fragmentation functions with the scale \i evolved up to 62 GeV using 
the Altarelli-Parisi equation (|2|). The average value of the momentum fraction for b — > B c is 
<z>= 0.68 at [L — 6.2 GeV, and it shifts downward to <z>= 0.61 when the scale is evolved 
up by a factor of 10. The process b — > B* has a slightly harder distribution with <z>= 0.73 
at jj, — 6.2 GeV and <z>= 0.66 at fi = 62 GeV. The fragmentation functions for c — ► B c 
and c — > B* are shown in Figure 3. The solid curves are the initial fragmentation functions 
at the scale \x = 11.2 GeV, where the running coupling constant has the value a s (fi) = 0.17. 
The fragmentation functions at the scale n = 112 GeV are shown as dotted curves in Figure 
3. The average value of the momentum fraction z is initially 0.51 for c — > B c and 0.55 for 
c — > B*. When the scale is evolved up by a factor of 10, they shift downward by about 9%. 
The evolution of the fragmentation functions would shift the momentum distributions for 
the B c and B* produced in Z° decay calculated by Chang and Chen [|] to slightly smaller 
values of z. 

The evolution equation (|2|) implies that at leading order in a s , the fragmentation 

probability Jq 1 dzD(z, /x) does not evolve with the scale /i. The evolution merely shifts the z 

distribution to smaller values of z. Therefore the fragmentation probabilities can be used to 

provide simple estimates of production rates. We need only multiply the production rates for 

bottom and charm quarks with energies much larger than their masses by the appropriate 

fragmentation probabilities. The fragmentation probabilities for the production of the B c 
are 

1 , n f s 2a s (/io) 2 |i?(0)| 2 / m c \ 

o dZ D "*&>*>) = 27^ f l^T^J ' (16) 
where the function f(r) is 

8 + 13r + 228r 2 - 212r 3 + 53r 4 r(l + 8r + r 2 - 6r 3 + 2r 4 ) , . , 

= + J (T^r - ] ^ r) ■ (17) 
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This function is almost constant, varying only from 0.496 at r = 0.23 to 0.464 at r = 0.77. 
The fragmentation probabilities for the production of the B* are 



\ n ( \ 2 Qs (/i ) 2 |/?(0)| 2 . 

/, w ,^« 5 m, d9) 

Jo ziTTin^ \mb + m c J 

where the function g(r) is 

24 + 109r - 126r 2 + 174r 3 + 89r 4 r(7 - 4r + 3r 2 + 10r 3 + 2r 4 ) , . . 

»w = ^(137^ + J l0E(r) - (20) 



It varies from 0.702 at r = 0.23 to 0.400 at r = 0.77. The fragmentation probability flT8"| ) 
agrees with the result for e + e~ — ► B*bc M obtained by Clavelli. The numerical values for 
the fragmentation probabilities are 2.2 x 10~ 4 for b — > B c , 3.1 x 10~ 4 for b — ► B*, 4.6 x 10 -6 
for c — > S c , and 3.9 x 10~ 6 for c — > 5*. The fragmentation probabilities for c are smaller 
than those for 6 by almost two orders of magnitude, because they are proportional to 1/ ml 
instead of 1/m 3 . 

It is interesting to compare our perturbative QCD fragmentation functions with the 
phenomenological Peterson fragmentation function |[T| which is often used to describe the 
fragmentation of a b quark into a hadron. For the case of a be hadron H, it is proportional 
to the energy denominator for the transition of a high momentum b into a be pair with zero 
relative velocity and a c: 

where e = m 2 /m 2 and iV# is a phenomenological parameter. The Peterson fragmentation 
function has the same qualitative behavior as the functions (|9|) and (|T2|) predicted by pertur- 
bative QCD. The main difference is that the QCD calculation predicts the normalizations of 
the fragmentation functions, while they are given by arbitrary phenomenological parameters 
Njj in the Peterson formula. The Peterson formula (^l|) also has the same shape for all 

bc hadrons H, while the precise shape of the QCD fragmentation functions depends on the 
orbital and spin quantum numbers. 
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The only inputs into our perturbative QCD calculations of the fragmentation func- 
tions are the QCD coupling constant a s , the quark masses and m c , and the nonrelativistic 
radial wavefunction at the origin R(0). The potential model values for R(0) should be quite 
reliable, since these models are tuned to reproduce the spectra of the charmonium and bot- 
tomonium systems and the be system represents an intermediate case. Potential models 
should also give accurate results for ratios of quark masses such as r = m c / (m& + m c ) . The 
largest uncertainty in the fragmentation functions for b — ► B c , B* should therefore come from 
the overall factor of 1/m 3 . Allowing for a variation of 0.2 GeV in the value m c = 1.5 GeV, 
the uncertainty in the normalization is an overall multiplicative factor of 1.5. Similarly, the 
largest uncertainty in the fragmentation functions for c — ► B c , B* comes from the overall fac- 
tor of 1/m 3 . Allowing for a variation of 0.3 GeV in the value mj = 4.9 GeV, the uncertainty 
is an overall multiplicative factor of 1.2 in the normalization. 

The fragmentation functions Di^ Bc (z, fi) and D c ^b c { z i aO can be used to calculate 
the rate for the direct production of B c in any high energy process. To obtain the total 
production rate for B c , we also need the direct production rates for all the other be states 
below the BD threshold, since they all decay ultimately into B c with branching fractions 
near 100%. We have calculated the fragmentation functions for the production of the 3 Si 
state B*. Our fragmentation functions can also be applied directly to the production of 
the first radial excitations of the B c and B* simply by changing the wavefunction at the 
origin. From potential model calculations ||10|| , the value for the first radial excitation of the 
S-wave states is |i?(0)| 2 = (0.99 GeV) 3 . Thus the production rates for these states should 
be about 60% of those for the B c and B*. The second radial excitation is probably above 
the BD threshold, so its branching fraction into B c will be negligible. Adding up all the 
S-wave contributions to the B c production rate, the total fragmentation probabilities are 
about 8.5 x 10~ 4 from the splitting of the b and 1.4 x 10 -5 from the splitting of the charm 
quark. Our lower bound of 8.5 x 10~ 4 for the probability for a high energy 6-quark to form 
a be bound state is consistent with the Monte Carlo estimates obtained in Ref. |J. 

An accurate calculation of the total production rate for B c could probably be obtained 
by including also the two sets of P-wave be states below BD threshold. In calculating the 
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production of the P-wave states, there are two distinct contributions that must be included 
at leading order in a s . The P-wave state can arise either from the production of a collinear 
be pair in a color-singlet P-wave state, or from the production of a collinear be pair in a 



color-octet S-wave state p2[ . There are two nonperturbative parameters that enter into the 
calculation. In addition to the derivative of the radial wavefunction at the origin for the 
P-wave state, we also require the probability density at the origin for the b and c to be in a 
color octet S-wave state accompanied by a soft gluon. 

Our results should be useful for assessing the prospects for the discovery of the B c in 
Z° decay at LEP or in pp collisions at the Tevatron. The branching fraction for Z° — > bb is 
approximately 15.2%. Multiplying by the total fragmentation probability for the production 
of S-wave be states, we obtain a lower bound on the inclusive branching fraction for the 
production of B c of 1.3 x 10~ 4 . The contribution from Z° — > cc followed by fragmentation 
of the charm quark or antiquark is smaller by about a factor of 80. At the Tevatron, the 
cross section for the production of charm quarks at large pt is a little larger than for bottom 
quarks, but not enough to make up for the difference in the fragmentation probabilities. 
Thus the dominant production mechanism for be states will again be b fragmentation. The 
cross section for the inclusive production of B c with px > 20 GeV at the Tevatron can be 
estimated by multiplying the cross section for the production of b quarks with px > 20 GeV, 
which is measured to be about 400 nb fI5| , by the fragmentation probability of 8.5 x 10~ 4 . 
Simple estimates of the production rates at future colliders can be obtained in the same way. 

We have pointed out in this paper that the dominant mechanism for the direct produc- 
tion of high energy be mesons is fragmentation, the production of a high energy b or c quark 
followed by its splitting into the the be state. We calculated the fragmentation functions 
D(z) for production of the S-wave states B c and B* to leading order in a s . These fragmen- 
tation functions are universal, so they can be used to calculate the rates and momentum 
distributions for the direct production of B c and B* in any high energy process. 

This work was supported in part by the U.S. Department of Energy, Division of High 
Energy Physics, under Grant DE-FG02-91-ER40684. 
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Figure Captions 



1. Feynman diagram for b* — > bcc which contributes to the fragmentation of b into B c 
and B*. 

2. Fragmentation functions D^ Bc (z, /i) and Di_> B * (z, ji) as a function of z at /i — 6.2 GeV 
(solid line) and fx — 62 GeV (dotted line). 

3. Fragmentation functions D c ^b c (z, fj) and D C ^ B * (z, /-*) as a function of z at /i = 11.2 GeV 
(solid line) and ji — 112 GeV (dotted line). 
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